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Preface 


There  are  14  figures  within  this  report,  many  with  a  level  of  complexity 
demanding  a  full-page  presentation.  To  avoid  the  disruption  to  the  text 
that  the  interspersing  of  multiple  full-page  illustrations  causes,  all  figures 
are  grouped  together  just  before  the  appendices.  It  is  hoped  that  the  reader 
does  not  find  this  to  be  an  inconvenience. 

Again,  with  regard  to  the  figures — all  surface  representations  were  produced 
by  using  the  analytic  surface  model  code  to  generate  an  array  of  surface- 
defining  points.  These  points  were  joined  by  the  plotting  routines,  for  con¬ 
stant  values  of  the  parametric  variables,  with  straight  line  segments.  Any 
perceived  surface  blemishes  are  attributable  to  the  plotting  routines  and  the 
simple  fashion  in  which  they  were  implemented  by  the  author. 

It  is  regretted  that  the  press  of  circumstances  has  not  permitted  the  de¬ 
velopment  of  any  truly  complicated  surfaces.  The  relative  simplicity  of  the 
modeled  shapes  should  not  be  taken  to  imply  any  undue  limitations  on  the 
capacity  of  the  model. 


Michael  J.  Vrabel 


There  are  more  ways  than  one  to  skin  a  cat. 


R.  A.  Habas.  Morals  for  Moderns 
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1  Analytic  Surface  Model  -  Theory 


1.1  Introduction 


The  ability  to  generate  high-fidelity  surface  models  is  fundamental  to  the 
solution  of  many  classes  of  problems  in  three  dimensions.  The  traditional 
approach  to  the  modeling  of  a  complex  surface  is  based  on  a  polynomial  fit 
to  an  array  of  surface-defining  points.  This  approach  can  suffer  from  sig¬ 
nificant  shortcomings.  In  particular,  depending  on  the  characteristics  of  the 
polynomial  fitting  algorithm,  the  exact  shape  of  the  resultant  surface  can 
be  difficult  to  predict  and  control.  More  specifically,  a  polynomial-generated 
surface  can  result  in  a  level  of  surface  irregularity  that,  while  not  always  visu¬ 
ally  apparent,  is  nevertheless  readily  detectable  when  a  surface  characteristic 
such  as  curvature  is  calculated. 


For  shapes  possessing  exact  analytic  descriptions,  these  shortcomings  do  not 
exist.  Exact  analytic  models  of  structures  are  generally  viewed  as  restricted 
to  objects  of  very  simple  shape.  To  be  meaningful,  such  an  approach  must 
extend  to  arbitrarily  complex  shapes.  This  report  develops  a  technique  for 
the  piece-wise  analytic  modeling  of  arbitrary  surfaces.  Succinctly,  the  ap¬ 
proach  is  to  design  a  surface  based  on  a  lofting  procedure  and  then  perform 
a  topological  transformation  on  the  lofted  model.  Lofting  refers  to  a  pro¬ 
cedure  for  creating  a  structure  by  “stretching”  a  surface  over  two  or  more 
parallel  structural  cross  sections.  Two  seminal  papers1,2  provided  the  stim¬ 
ulus  for  the  development  of  this  model— the  first  by  developing  the  basis  for 
an  analytic  lofted  surface,  and  the  second  by  demonstrating  the  flexibility 
of  the  superquadrics  for  shape  modeling.  It  should  be  made  clear  from  the 
onset  that,  while  the  title  of  this  report  indicates  a  model  restricted  to  aero¬ 
dynamic  shapes,  the  model  is  general  and  can  be  used  for  a  wide  range  of 
shapes.  It  is  for  aerodynamic  shapes  with  their  smooth  (and  smooth  is  used 
here  in  its  colloquial  sense)  transitioning  surfaces  that  the  model  is  most 
easily  adapted.  Included  within  the  report  is  the  full  theoretical  framework, 


:R.  L.  Barger  and  M.  S.  Adams,  “Semianalytic  Modeling  of  Aerodynamic  Shapes,”  National  Aeronautics 
and  Space  Administration,  NASA  TP-2413,  1985. 

A.  H.  Barr,  “Superquadrics  and  Angle-Preserving  Transformations,”  IEEE  CO, A,-  A  Vol  1  No  1  Jan 
1981.  ’  '  ’ 
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a  users  guide  for  a  FORTRAN  implementation  of  the  model,  a  source  code 
listing,  and  numerous  examples. 


1.2  The  Barger  Model 

The  Barger  and  Adams  paper  describes  a  lofting  technique  for  the  analytic 
modeling  of  classes  of  smoothly  varying  surfaces.  A  surface  is  defined  as  a 
sequence  of  subsurfaces  formed  by  the  transitioning  between  bounding  cross- 
sectional  shapes.  The  equation  for  a  subsurface  in  Cartesian  coordinates  is 

y{x,t)  =  S(x){[l  -  T(x)]yi(t)  +  T(x)y2(t)}  (l) 

z(x,t)  =  S(x){[l  -  T(x)]zi(t)  +  T(x)z2{t)}  (2) 

where  yi(t),  zi(t),  y2(t),  and  z2(t)  are  analytic  functions  in  parametric  form 
of  the  bounding  cross  sections,  T(x)  is  the  transition  function  governing  the 
shape  of  the  subsurface  and  varying  from  1  to  0  as  x  varies  from  0  to  1, 
5(x)  is  the  size  function  providing  an  independent  function  scaling,  t  is  a 
parametric  variable,  and  x  is  the  normalized  (over  a  subsurface)  longitudinal 
coordinate.  The  condition  imposed  on  the  transition  function  at  the  interface 
between  two  subsurfaces  (A  and  B)  where  continuity  is  required  is 

T%  (1)  =  Tq  (0)  =  0  (3) 

where  N  defines  the  order  of  the  derivative  for  which  continuity  is  to  be 
enforced.  The  further  conditions  on  T(x )  are  that  it  be  single- valued  and 
monotonic.  Since  a  surface  describable  by  an  arbitrary  transition  function 
can  be  divided  into  elements  for  which  monotonicity  can  be  enforced,  this 
latter  condition  poses  no  undue  surface  modeling  limitation. 

Although  equations  (l)  and  (2)  indicate  that  the  size  function  £(z)  operates 
over  a  single  subsurface,  where  multiple  continuous  elements  are  defined, 
5(x)  will  operate  over  each  continuous  subset.  In  this  way,  interelement 
continuity  is  not  affected.  The  conditions  on  5(x)  are  that  it  be  smooth, 
single- valued,  and  non-negative. 


1.3  Analytic  Term  Definition 

The  previous  section  provides  a  framework  for  defining  a  class  of  analytic 
surfaces.  It  is  appropriate  at  this  point  to  both  generalize  terms  within  the 
lofted  surface  equations  and  change  notation,  substituting  for  t  and  x  the 
parametric  variables  u  and  v.  This  latter  change  permits  the  three  compo¬ 
nents  of  the  surface-defining  vector  to  be  expressed  as  dependent  variables: 

yi{u,v)  =  5„(u){[l  -  Tv(u)]yi(u)  +  Tv(y)y2(u)}  (4) 
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(5) 

(6) 


z,(u,v)  =  St{v){[l  -  T,{ v)}Zl{u)  +  Tz{v)z2{u)} 

Xi(v)  -=  V 

The  reasons  for  these  changes  will  become  apparent. 

To  be  a  viable  model,  for  the  generalized  transition,  size,  and  cross-section 
functions,  expressions  must  be  substituted  that  are  a  function  of  and  can 
be  readily  correlated  with  geometric  shape.  Rather  than  requiring  the  user 
to  define  unique  expressions  for  each  geometric  construct,  this  section  in¬ 
tends  to  demonstrate  a  set  of  equations  general  enough  to  be  used  without 
modification  for  defining  a  broad  range  of  shapes. 

1.3.1  Transition  Function 

A  superellipse  is  defined  as 


y  =  A\  sin6  u 
z  =  A2  cos6  u 

where  0  <  u  <  2n  and  e  >  0.  To  generate  a  four  quadrant  plot 
these  equations  must  be  modified: 


(7) 

(8) 

(see  fig.  1), 


y  =  Ai  sin  u  |  sin  u  |6  1  (9) 

Z  —  A2COSU  I  COSK  |c_1  (10) 

To  preserve  the  functional  simplicity  of  the  original  definition  and  allow  its 
use  in  all  subsequent  operations,  the  following  rule  is  adopted:  the  sign  of 
the  sin  and  cos  terms  is  to  be  taken  out  from  under  the  exponent.  This  rule 
makes  equations  (7)  and  (9)  and  equations  (8)  and  (10)  equivalent  and  can 
be  carried  through  to  define  a  proper  set  of  function  derivatives. 

Equations  (7)  and  (8)  can  be  combined: 


y  =  a,  [1  -  (z/a2)2/*]c/2 

(11) 

If  for  e  one  substitutes 

e  =  a!  +  a2z°3  +  a4zab 

(12) 

where  alt  ^5  are  constants,  then  equation  11  can 
equation  that  can  be  used  as  a  transition  function: 

be  transformed  into  an 

o* 

0 

1 

rH 

II 

& 

(13) 
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where 


(14) 


<3  =  2/  (fl]  +  tt2Va3  +  <Z4l>a'J) 

An  alternate  form  for  Q  is 

Q  =  2/  [ai  +  a2ua3  +  a4sina:,(7ru)]  (15) 

This  expression  generates  a  large  family  of  monotonic  curves  meeting  the 
requirements  that  T(0)  =  1  and  T(l)  =  0,  and  the  conditions  imposed  by 
equation  (3)  to,  at  least,  the  second  derivative  level  at  both  v  =  0  and 
v  =  1.  This  is  demonstrated,  for  a  select  sample,  in  figures  2  through  4.  The 
transition  functions  demonstrated  in  these  figures  are  intended  to  meet  the 
indicated  boundary  value  conditions.  Where  less  rigorous  conditions  are  to 
be  met,  quite  different  curve  shapes  are  possible.  To  meet  first-derivative 
continuity  conditions  the  limits  on  c  are  e  <  1.4  about  v  =  0  and  e  >  2.9 
about  v  —  1.  To  meet  second-derivative  conditions,  e  <  0.8  about  v  =  0 
and  e  >  4.9  about  v  =  1.  There  is  little  loss  in  the  generality  of  the  lofted 
surface  model  if  equation  (13)  is  made  the  transition  function.  The  Barger 
model  includes  the  same  transition  function  for  the  y  and  z  components  of 
the  vector  defining  the  lofted  surface.  A  useful  generalization  is  to  allow  the 
transition  function  to  be  different  for  y  and  2.  This  is  indicated  in  equations 
(4)  and  (5)  and  is  included  in  the  codes  listed  in  appendix  A  and  the  final 
equation  summary. 

1.3.2  Cross  Section 

Figure  1  demonstrates  the  diversity  of  shapes  possible  with  the  simple  func¬ 
tion  defining  the  superellipse.  This  poses  the  possibility  that  some  variation 
on  the  superellipse  may  prove  useful  in  defining  a  generalized  equation  for 
modeling  cross  sections.  Central  to  this  task  is  the  recognition  that  the  rela¬ 
tionship  between  the  parametric  variable  and  the  y-z  location  of  the  function 
is  well  defined  for  four  values  of  u  independent  of  e: 


u  = 

0 

:  +zaxis 

U  —  7r/2 

:  +yaxis 

u  = 

7T 

:  — zaxis 

u  —  37t/2 

:  — yaxis 

Consequently,  a  more  general  (and  more  useful)  expression  can  be  developed: 

y  =  Ai  sin6  (mu  +  uy)  (16) 

2  =  A2  cos6  (mj  +  uz)  (17) 
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where  m,  n ,  uy,  and  uz  are  constants  that  can  be  redefined  whenever  the 
following  relationships  are  satisfied: 

mu  +  uy  =  7tkm/2 

nu  +  uz  =  nkn/2 

for  km,kn  =  0,1,2,....  The  resultant  function  will  be,  at  least,  analytic 
in  each  interval.  Further  flexibility  results  from  making  e  a  function  of  the 
parametric  variable: 

c  =  61  +  biqh:i  (18) 

q  =|  sin  (ju  +  u9)  | 

v 

where  the  constants  associated  with  the  argument  of  q  are  governed  as  in 
equations  (16)  and  (17). 

These  equations  define  curves  with  fixed  locations  with  respect  to  the  co¬ 
ordinate  system  origin  and  axes.  To  allow  for  rotation  and  translation  in 
a  plane,  an  affine  transformation  is  performed.  The  resultant  cross-section 
equations  are 


y  -  Az  { [Aj  sine  (mu  +  uv)  +  ytl 

cos  xp  + 

A2  cos6  (nu  +  uz )  +  ztl 

sin  +  y,. 

r  r  f 

(19) 

z  =  A4  (  -  [Ai  sin  (mu  +  uy)  +  yix 

sin  xp  + 

A2  cos6  (nu  +  uz)  +  ztl 

cos  ip j  +  Zt2 

(20) 

where  ytl,  yt2,  ztl,  and  zt2  define  the  components  of  two  independent  dis¬ 
placement  vectors  and  xp  defines  the  angle  of  rotation  about  the  origin.  The 
translation  and  scaling  redundancies  included  in  these  equations,  by  allow¬ 
ing  the  corresponding  operations  both  before  and  after  a  rotation,  facilitate 
model  use. 

The  rotation  and  translation  terms,  combined  with  the  parameter  redefini¬ 
tion  capability,  permit  a  cross-section  model  fragmented  into  an  unlimited 
number  of  open  and  closed  curves.  Such  a  capability,  combined  with  the 
lofted  surface  model,  allows  the  definition  of  multibranched  bodies. 

1.3.3  Size  Function 

The  size  function  is  intended  to  provide  an  independent  scaling  of  the  lofted 
surface.  Although  the  size  function  is,  in  a  sense,  redundant,  its  presence 
provides  a  measure  of  model  flexibility  difficult  to  duplicate  with  just  the 
cross-section  and  transition  functions.  Two  useful  forms  for  the  size  function 
are 

m 

‘S’H  =  J2  ClnVC2n  (21) 

n=  1 
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and  an  alternative, 


S(v) 


E  (ClnU  +  C2n)C3' 

.n=l 


C4 


(22) 


where  m  is  generally  less  than  or  equal  to  4,  0  <  v  <  1  ,  and  Cini...i3n,C4  can 
take  on  any  value  consistent  with  the  restrictions  imposed  upon  S(v). 


As  with  the  transition  function,  allowing  for  independent  size  functions  for 
the  y  and  z  surface-defining  vector  components  will  prove  useful. 


1.4  Topological  Transformation 


The  strength  of  a  lofted  model,  namely  its  ability  to  define  a  shape  through 
the  specification  of  a  limited  number  of  cross  sections,  is  also  its  weakness. 
At  the  model’s  ends,  that  is  at  v  =  0  and  v  =  1,  shape  definition  is  con¬ 
strained  by  the  requirement  that  the  end-defining  cross  sections  be  planar 
and  normal  to  the  model  axis  (x  axis).  This  limitation  can  be  circumvented 
by  performing  a  topological  transform  on  the  lofted  model.  This  transfor¬ 
mation  is  to  operate  upon  the  x  component  of  the  surface-defining  vector, 
substituting  for  the  simple  expression  of  equation  (6)  one  that  is  a  function  of 
both  u  and  v.  Properly  implemented,  such  a  transformation  allows  complete 
freedom  in  specifying  an  out-of-plane  shape  for  any  model  cross  section. 


Because  of  the  success  of  the  lofted  model,  the  topological  transformation 
will  employ  similar  equations.  For  equation  (6),  one  substitutes 

x(u,v)  —  xt(v)  +  xt(u,v)  (23) 


where  xt(u,v)  is  yet  to  be  defined.  Equations  (4)  and  (5)  are  left  unchanged: 

y{u,v)  =  yi(u,v)  (24) 

z(u,v)  =  zi(u,v)  (25) 

Analogous  to  the  basic  lofted-model  equation,  xt(u,v)  is  defined  as 

xt{u,  v)  =  (l  -  x^(u))  Gf(u,  v )  +  xfc(v)Gb(u,  v )  (26) 

where  xin(v )  is  xi[v)  normalized  over  an  appropriate  interval  (either  the  in¬ 
terval  between  adjacent  cross  sections  or  the  length  of  the  modeled  structure) 
and  h  is  a  constant  that  defines  the  rate  at  which  xt(u,v)  transitions  between 
the  two  terms  of  this  equation.  The  terms  Gf(u,v)  and  Gt(u,v)  are  defined 


as 

Gf{u,v) 


h 


n 


h 


f  2 


Vfi  l2/c' 


le//2 


+  h 


/3 


1-  \hf4-  Zfi  I2/6/ 


le//2 
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+  fi-yio 
(27) 


Gb(u,v)  =  hbl  [l-  |  hb2  -  ybl  |2/et]e‘/2  +  hb 3  [l-  |  hH  -  zbl  |2/e,]et/2  +  hbw 

(28) 

0  <  yfM  <  1  (29) 

0  <  zfM  <  1  (30) 

with 

ef  =  ^/5  +  hf6yj{7  +  hf8zhji0  (31) 

*6  =  hb  5  +  hbsyft7  +  /i6826fc/°  (32) 

where  yfibl  and  z/,«  are  the  normalized  y  and  z  coordinates  of  cross-sections  / 
and  b  front  and  back  and  hfib i.^io  are  constants.  The  similarity  between 
equations  (26)  through  (32)  and  equations  (4),  (5),  (12),  and  (13)  should  be 
noted. 


To  increase  the  flexibility  of  equation  (26),  parameter  switching  can  be  in¬ 
cluded.  Since  section  1.3.2  permits  the  definition  of  each  cross  section  in 
terms  of  a  composite  of  discrete  curve  segments,  equation  (26)  can  em¬ 
ploy  parameter  switching  for  each  discrete  segment  or  for  clusters  of  these 
segments. 


1.5  Three-Dimensional  Model 

Equations  (23),  (24),  and  (25)  are  the  components  of  a  surface-defining  vec¬ 
tor.  This  vector  contains  69  user-defined  constants.  With  parameter  switch- 
in8>  the  total  number  of  constants  is  effectively  limited  only  by  the  structure 
of  the  model-implementing  computer  code.  This  number  presently  exceeds 
15,000.  The  obvious  implication  is  that  the  complexity  of  the  surface  capable 
of  being  modeled  is  limited  by  the  user’s  capacity  to  define  the  appropriate 
parameters. 

A  compact  final  summary  of  the  model  equations  follows. 
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yi{u,v)  =  Sy(v)  {(1  -  Tf(t>)]  yi(u)  +  Ty(v)y2{u)}  (33) 

zi(u,v)  =  Sz(v)  {[1  -  Tz(v)}  zi(u)  +  Tz(v)z2(u)}  (34) 

xi(v)  =  v  (35) 

TyAv)  =  ~  v9)'IQ  (36) 

<3  =  2/  (a,,.,  +  a2t,^as».*  +  (37) 

Sy,z(v)  =  ciyiZ  +  C2y,zvCzv'z  +  ciytZvChy'z  (38) 


yi  2(u)  =  A3  |  [Ai  sine  (mu  +  uy)  +  y*,]  cos  ip  +  [.A2  cose  ( nu  +  uz)  +  ztl j  sin  V’}  +  Vt2 

(39) 

z\  2  (u)  =  A4  |  —  \Ai  sine  (mu  +  uy)  +  y(l]  sin  ip  +  [a2  cose  (nu  +  uz )  +  ztl  j  cos  ip  j  +  zt2 

(40) 


e  —  61  +  629^  (41) 

q  =|  sin  (ju  +  uq)  |  (42) 

y(u,v)  -  yi(u,v)  (43) 

z(u,v)  —  zi(u,v )  (44) 

x(u,v)  =  x;(u)  +  ^1  -  xfc'j  Gf(u,v )  +  x^Gb(u,v)  (45) 


r  *16^/2  r  “j  £//2 

Gf(u,  v)  =  hfi  1—  |  hj2  —  yfi  \2/et  +  hf3  1 1-  |  hf4  -  zfi  |2/e'  j  +  hfl0 

(46) 

Gb(u,v)  —  hbi  [l-  |  hb2  -  ybt  |2/e»]  ^  +  hb3  [l-  |  hb4  -  zbt  |2/e*>]  +  hbw 

(47) 

0  <  xin,yftbi,zfibi  <  1  (48) 

e/,6  =  ^/,65  +  hfibey/bi7  +  hfibgz^bi°  (49) 


1.6  Conclusion 

A  piece-wise  analytic  function  capable  of  defining  a  wide  range  of  geometric 
surfaces  can  be  formed  by  combining  a  modified  superellipse  with  a  topo¬ 
logically  transformed  lofted  model.  The  final  form  of  the  model’s  equations 
represents  a  reasonable  compromise  between  the  desire  to  keep  the  model 
simple  (and  hence  easy  to  use)  and  to  ensure  the  capability  of  modeling 
complicated  surfaces.  The  interface  with  this  function  is  an  array  of  user- 
defined  constants.  Because  these  parameters  are  easily  correlated  with  the 
shape  of  the  modeled  structure,  shape  definition  and  alteration  is  readily  im¬ 
plemented.  For  applications  restricted  to  the  well-behaved  surfaces  of  exact 
analytic  models,  this  function  opens  a  host  of  new  geometries. 
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2  Analytic  Surface  Model  -  Code 


2.1  Code  User’s  Manual 

The  user  s  manual  is  a  guide  to  the  use  of  the  model-implementing  computer 
code  (FORTRAN  source  module  listed  in  app  A).  The  code  is  designed  to 
make  optimum  use  of  the  capabilities  implicit  in  the  theoretical  development 
of  the  previous  sections.  This  means  that  more  flexibility  is  built  into  the 
code  than  is  explicitly  indicated  by  the  theoretical  treatment.  This  will 
become  apparent  within  the  body  of  the  user’s  manual. 

All  input  to  the  code  is  via  a  namelist  statement.  Each  term  within  the 
namelist  statement  is  defined  below. 


2.1.1  xl 


xl  relates  to  the  longitudinal  (x-axis)  length  of  the  lofted  model.  The  code 
has  the  capability  of  treating  a  maximum  of  10  independent  structures  in 
a  single  pass.  These  structures  can  be  joined  at  common  cross  sections 
to  form  a  composite  body.  There  are  several  advantages  to  using  multiple 
independent  structures  to  define  a  single  body: 

•  Each  independent  structure  can  be  developed  and  tested  separately, 
reducing  the  modeling  burden  on  the  user. 

•  Structures  that  are  not  monotonic  in  x  can  be  modeled. 

•  The  parameters  within  the  size  and  transition  functions  can  be  changed 
over  the  length  of  the  composite  structure. 

xl  is  the  x-axis  location  of  all  model  cross  sections  (station  cuts)  listed  in 
order,  xl  must  be  monotonic  in  x  within  an  independent  structure  but  need 
not  be  monotonic  between  structures.  Where  independent  structures  are 
to  be  joined  at  a  common  cross  section,  the  x-axis  location  of  that  station 
cut  must  appear  twice  in  the  xl  sequence.  When  using  a  composite  body, 
one  must  exercise  caution  to  ensure  that  continuity  conditions  are  met  at  all 
interfaces. 
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2.1.2  ncount 


ncount  is  the  number  of  independent  structures. 


2.1.3  ns 

ns  is  the  number  of  model  segments  defined  by  station  cuts  within  each 
independent  structure  (ns  =  number  of  station  cuts  less  1). 

2.1.4  anx 

anx  relates  to  the  density  of  calculated  surface  points  and  is  the  surface  grid 
count  in  the  x  direction  for  each  independent  structure.  For  the  x  direction, 
surface  grid  density  is  constant  within  individual  independent  structures. 


2.1.5  mid 

mid  relates  to  the  relative  placement  of  the  calculated  surface  point  for  the 
x  coordinate  of  the  position  vector.  For  mid  =  0,  the  first  point  falls  on  the 
first  grid  line,  and  subsequent  points  on  subsequent  grid  lines.  The  first  grid 
line  is  coincident  with  the  corresponding  first  station  cut.  For  mid  ^  0,  each 
point  falls  midway  between  grid  lines. 

2.1.6  npoint 

npoint  is  the  surface  grid  count  transverse  to  the  x  direction  for  each  inde¬ 
pendent  block. 


2.1.7  ang 

ang  determines  the  relative  point  distribution  about  the  perimeter  of  each 
station  cut.  Each  station  cut  can  be  composed  of  a  maximum  of  15  curve 
segments.  The  summation  of  ang  for  all  15  segments  (or  utilized  fraction 
thereof)  must  equal  360.  The  grid  spacing  is,  generally,  not  uniform  within 
each  curve  segment  but  tends  to  cluster  more  densely  in  regions  where  cur¬ 
vature  is  changing  most  rapidly. 


2.1.8  nseg 

nseg  is  the  number  of  curve  segments  contained  within  each  station  cut  for 
a  maximum  of  25  station  cuts. 
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2.1.9  uy 


uy  is  the  starting  value  in  degrees  of  the  parametric  variable,  u,  for  the  y  co¬ 
ordinate  of  each  individual  segment  composing  a  station  cut.  It  corresponds 
to  uy  in  equation  (39). 

2.1.10  uz 

uz  is  the  corresponding  value  for  the  z  coordinate,  and  for  the  present  code 
structure  must  be  set  identically  equal  to  uy.  It  corresponds  to  uz  in  equation 
(39). 

2.1.11  alim 

alim  is  the  corresponding  stopping  value  in  degrees  of  the  parametric  vari¬ 
able,  u,  for  both  y  and  z  coordinates  independent  of  the  values  of  am  and 
an.  alim  can  be  made  larger  or  smaller  than  uy  and  uz.  For  values  of  alim 
or  uy  and  uz  spanning  360  degrees,  use  an  appropriate  negative  value  for  the 
stop  or  start  point  so  that  the  absolute  value  of  alim  minus  uy  remains  the 
angular  extent  of  the  spanned  interval.  To  calculate  the  appropriate  value 
of  alim,  temporarily  treat  parameters  am  and  an  as  equal  to  1.  In  this 
way  the  full  angular  extent  spanned  in  the  definition  of  a  cross  section  or 
cross-section  segment  can  be  divorced  from  the  0  to  360  degree  span  of  the 
parametric  variable,  u. 

2.1.12  uq 

uq  corresponds  to  uq  in  equation  (42)  and  is  analogous  to  uy  and  uz.  The 
range  spanned  by  the  argument  of  the  sin  function  in  equation  (42)  is  gov¬ 
erned  by  alim. 

2.1.13  am 

am  corresponds  to  m  in  equation  (39). 

2.1.14  an 

an  corresponds  to  n  in  equation  (39). 

2.1.15  aj 

aj  corresponds  to  j  in  equation  (42). 
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2.1.16 


bl  b2  b3 


2.1.17 

2.1.18 

2.1.19 

2.1.20 

2.1.21 

2.1.22 

2.1.23 

2.1.24 

2.1.25 


bl,  b2,  and  b3  correspond  to  61,62,  and  63  in  equation  (41). 

a  b 

a  and  b  correspond  to  Aj  and  A2,  respectively,  in  equation  (39). 

aa  bb 

aa  and  bb  correspond  to  A3  and  ^44,  respectively,  in  equations  (39)  and  (40). 

psi 

psi  corresponds  to  0  in  equation  (39). 


zt  yt 

zt  and  yt  correspond  to  zt2  and  y<2,  respectively,  in  equations  (39)  and  (40). 

ztp  ytp 

ztp  and  ytp  correspond  to  ztl  and  t/f, ,  respectively,  in  equation  (39). 

al  a2  a3  a4  a5 

al,. ..,  a5  correspond  to  alv  ..,a5,  respectively,  in  equation  (37),  where  the 
transition  function  is  identical  for  both  y  and  z. 

ayl  ay  2  ay  3  ay  4  ay  5 

ayl,. . .,  ay5  correspond  to  ai„, . . . ,  a^y  in  equation  (37),  where  the  transition 
function  is  different  for  y  and  z. 

azl  az2  az3  az4  az5 

azl,. . .,  az5  correspond  to  a\z, . . .  ,a^z  in  equation(37),  where  the  transition 
function  is  different  for  y  and  z. 

cl  c2  c3  c4  c5 

cl,. . .,  c5  correspond  to  cj, . . .  ,c5  in  equation  (38),  where  the  size  function 
is  identical  for  both  y  and  z. 
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2.1.26 

2.1.27 

2.1.28 

2.1.29 

2.1.30 

2.1.31 

2.1.32 


cyl  cy2  cy3  cy4  cy5 

cyl,. . .,  cy5  correspond  to  cjy, . . . ,  cby  in  equation  (38),  where  the  transition 
function  is  different  for  y  and  z. 

czl  cz2  cz3  cz4  cz5 

czl,. . cz5  correspond  to  c\z , . . .  ,cbz  in  equation  (38),  where  the  transition 
function  is  different  for  y  and  z. 

h  hfl  hf2  hf3  hf4  hf5  hf6  hf7  hf8  hf9  hflO 

h,  hfl,...,  hflO  correspond  to  h,  hfu . . .  ,hflQ  in  equations  (46)  and  (49). 

hbl  hb2  hb3  hb4  hb5  hb6  hb7  hb8  hb9  hblO 

hbl,. . hblO  correspond  to  hbu  . .  ,,hbl0  in  equations  (47)  and  (49). 


ntop 

ntop(n)  =  1  is  used  to  show  that  a  topological  transform  is  to  be  performed 
upon  independent  structure,  n.  Otherwise,  ntop(n)  is  set  to  0. 

nwi 

nwi  determines  the  x  interval  over  which  the  topological  transform  operates 
within  an  independent  structure.  If  nwi  is  set  to  1,  then  the  topological 
transform  operates  over  the  entire  independent  structure.  For  nwi  set  to 
0,  it  operates  over  individual  segments  (segment  count  determined  by  ns) 
within  the  independent  structures. 

ngroup 

ngroup  defines  clustering  of  cross-section  curve  segments  for  use  by  topolog¬ 
ical  transform,  ngroup  is  composed  of  groups  of  five  numbers.  The  first  and 
third  are  the  cross-section  numbers  (same  value).  The  second  is  the  starting 
cross-section  curve  number.  The  fourth  is  the  stopping  cross-section  curve 
number  (the  second  and  fourth  are  permitted  to  be  the  same).  The  fifth  is 
a  requisite  0.  Within  each  five-unit  cluster  defined  by  ngroup,  all  previ¬ 
ously  defined  h  parameters  are  held  constant.  List  all  ngroup  parameters 
in  ascending  order. 
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2.1.33  list3  nblock 


list3  determines,  in  part,  the  format  of  the  code  output.  For  list3  equal  to 
0,  the  x  component  of  the  surface-defining  vector  is  stored  on  FORTRAN 
logical  unit  jx,  y  is  stored  on  jx+1,  and  z  is  stored  on  jx+2.  For  list3  not 
equal  to  0,  x,  y,  and  z  are  stored  on  logical  unit  jx  (3el5.6). 

nblock  allows  the  output  to  be  stored  on  separate  groups  of  FORTRAN 
logical  units  for  different  groups  of  independent  structures.  With  the  inde¬ 
pendent  structures  numbered  sequentially  and  in  ascending  order,  nblock 
is  the  independent  structure  number  of  the  last  contiguous  member  of  each 
independent  block.  Each  independent  block  is  incremented  by  3  with  jx  ini¬ 
tially  set  to  60  (see  list3).  This  arrangement  permits  any  combination  of 
nblock  and  list3. 


2.2  Model  Demonstration  —  Examples 

This  section  provides  some  guidelines  for  the  use  of  the  model.  These  guide¬ 
lines  primarily  demonstrate  the  characteristics  of  the  model  equations. 

2.2.1  Transition  Function 

In  the  theory  section,  the  discussion  of  the  requirements  imposed  on  the 
transition  function  for  interelement  continuity  (eq.  (3))  was  based  on  the 
most  general  interpretation  of  the  Barger  model.  Unfortunately,  the  require¬ 
ments  of  this  equation  can  lead  to  undesirable  structural  shapes  or  severe 
limitations  on  the  shape  of  the  modeled  structure.  For  this  reason,  examin¬ 
ing  the  basis  for  equation  (3)  is  worthwhile.  Consider  a  structure  formed  by 
two  elements:  a  (defined  by  cross  sections  1  and  3)  and  b  (defined  by  cross 
sections  3  and  2  ).  (Refer  to  the  original  NASA  paper  [ref  l]  for  details  of  the 
derivation  of  the  following  equation.)  To  ensure  first-derivative  continuity  at 
common  cross  section  3,  the  following  relationship  must  be  satisfied: 

y2  -  yz  dTb  _  t/3  ~  Vi  dTg 

X2  «X3  X3  Xj  dzE/j 

where  an  identical  expression  exists  for  the  z  component.  For  certain  values 
of  the  x  and  y  parameters,  this  relationship  can  be  satisfied  for  conditions 
other  than  dTb/dxb  =  dTa/dxa  —  0. 

•  If  y1  =  y2  =  y3'5  this  relationship  is  satisfied  for  any  form  of  7)  and 
Ta.  This  seemingly  trivial  example  is,  in  fact,  important  because  many 
situations  can  be  encountered  where  continuity  is  required  over  only 
a  part  of  a  surface  defined  by  three  (or  more)  identical  partial  cross 
sections. 
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•  If  (j/2  —  2/3) / (^2  ~  £3)  =  (3/3  —  Vi)/{^3  ~  %i),  then  the  relationship  of 
equation  (50)  is  satisfied  for  dT^/dxi,  =  dTa/dxa.  This  situation  is 
encountered,  for  example,  for  a  frustum  of  a  cone  where  the  common 
cross  section  is  located  anywhere  between  the  parallel  end  surfaces.  In 
general,  this  condition  is  valid  for  surfaces  of  constant  slope — a  much 
encountered  situation. 

The  same  relationships  can  be  shown  to  hold  true  for  higher  order  deriva¬ 
tives.  This  broadening  of  interelement  continuity  requirements  significantly 
increases  the  utility  of  the  analytic  model. 

The  form  of  the  exponent  in  equation  (37)  was  selected  because  of  its  ability 
to  transition,  in  both  a  complicated  and  readily  predictable  fashion,  among 
the  first  quadrant  curves  of  the  superellipse  as  v  goes  from  0  to  1.  In  lieu 
of  generating  reams  of  transition  function  curves  for  various  combinations 
of  exponent  parameters,  the  most  useful  tool  for  designing  an  appropriate 
transition  function  becomes  merely  a  family  of  curves  detailing  the  behavior 
of  for  a  broad  range  of  values  of  N.  This  is  done  in  figure  5. 


2.2.2  Cross-Section  Function 

The  best  way  to  demonstrate  the  performance  of  the  cross-section  function  is 
to  start  with  the  simplest  form  of  the  function,  demonstrate  its  performance, 
and  repeat  the  process  as  additional,  complicating  terms  are  included. 

Step  1 


y  =  M  sin6  (mu  +  uy )  (51) 

z  =  A2  cos  e(nu  +  uz)  (52) 


For  step  1,  A\  =  A2  =  1,  and  e  is  held  constant.  The  constants  uy  and  uz 
define  phase  terms  for  y  and  2.  The  consequences  of  setting  uy  ±  uz  are,  for 
the  more  complex  cross-section  models,  difficult  to  predict  and,  for  a  model 
based  on  intuitive  parameters,  therefore  difficult  to  exploit.  For  this  reason, 
this  inequality  is  not  permitted  by  the  code.  The  consequences  of  setting 
uv  7^  uz  in  equations  (51)  and  (52)  are  demonstrated  in  figure  6. 

The  constants  m  and  n  define  the  range  of  the  arguments  of  the  sin  and  cos 
functions  and,  unlike  with  uy  and  uz,  having  m  ^  n  can  prove  useful.  The 
consequences  of  setting  m  ^  n  are  also  demonstrated  in  figure  6,  along  with 
the  results  of  simultaneously  varying  both  sets  of  parameters.  At  this  point 
a  few  observations  are  appropriate. 
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•  The  generated  shapes  are  complicated. 

•  The  generated  shapes  are  not  particularly  interesting. 

•  Although  the  curve  shapes  should  be  predictable,  they  are  not  easily 
predictable. 

Step  2 


y  =  A\  sin6  (mu  +  uy )  (53) 

z  —  Ai  cos6  (nu  +  uz )  (54) 

e  =  &i  +  (5^) 

q  =|  sm  (ju  +  uq)  |  (56) 


In  this  step,  e  is  made  a  function  of  the  parametric  variable.  For  the  proper 
operation  of  equation  (55),  it  is  highly  desirable  to  constrain  q  to  the  range 
0  to  1;  hence  the  form  of  equation  (56).  The  similarity  in  form  between 
equation  (56)  and  equations  (53)  and  (54)  is  not  coincidental,  as  the  same 
considerations  apply.  Additionally,  the  similarity  between  equations  (55) 
and  (37)  is  intentional  since,  again,  the  same  considerations  are  involved. 
The  variable  e  functions  by  transitioning  among  the  family  of  curves  defined 
by  the  superellipse  as  the  parametric  variable  is  exercised  over  its  range. 
Examples  of  step  2  equations  are  given  in  figure  7.  To  constrain  the  number 
of  curve  shapes,  all  results  represent  variations  on  a  circle.  What  this  means 
should  be  obvious  from  the  curves. 

Step  3 

This  step  operates  off  the  same  set  of  equations  as  the  previous  step,  except 
all  parameters  are  allowed  to  change  whenever  the  following  conditions  are 
satisfied: 


mu  +  uy  —  nkm/2 

(57) 

nu  +  uz  =  Trkn/2 

(58) 

ju  +  uq  =  7tkj/2 

(59) 

where  km,kn,kj  =  0,1,2,....  Since  examples  of  the  model  at  this  step  can 
be  constructed  by  splicing  together  segments  of  the  curves  in  the  previous 
step,  no  specific  examples  are  given. 
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Step  4 

The  final  step  allows  arbitrary  rotation  and  translation  of  model  segments: 

y  =  A 3  1 1  A\  sine  (mu  +  uy)  +  ytl j  cos  ip  +  [a2  cose  (nu  +  uz )  +  ztl j  sin  ^  J  +  t h2 

(60) 

2  =  A4  {  -  [Ai  sine  (mu  +  uy)  +  yfl]  sin  ip  +  [a2  cose  (nu  +  uz)  +  ztl ]  cos  0}  +  zti 

(61) 

Examples  of  the  capability  of  the  final  cross-section  model  are  given  in 
figure  8. 

2.2.3  Size  Function 

The  size  function  provides  an  independent  scaling  of  the  lofted  surface.  By 
allowing  different  size  functions  for  y  and  z,  a  more  flexible  model  results. 
To  demonstrate  this,  a  right  circular  cylinder  is  to  be  operated  upon  by  a 
range  of  y  and  z  independent  size  functions.  The  results  are  demonstrated 
in  figure  9. 

2.2.4  Topological  Transformation 

The  topological  transformation  derives  its  name  from  the  manner  in  which 
it  operates  upon  the  lofted  model.  The  transformation  performs  an  out-of¬ 
plane  distortion  of  selected  model  cross  sections  with  the  consequent  alter¬ 
ation  in  surface  shape.  The  similarity  in  form  between  the  equations  of  the 
lofting  procedure  and  the  topological  transformation  precludes  the  need  for 
any  extended  discussion  of  model  behavior  except  to  note  that  parameters 
hf,b 2  and  hJM  are  generally  restricted  to  values  of  1  or  0.  An  example  of  a 
topologically  transformed  structure  is  given  in  figure  10.  The  starting  lofted 
model  was  two  identical,  joined,  right-angle  rectangular  tubes. 

2.2.5  Three-Dimensional  Model 

To  demonstrate  the  performance  of  the  analytic  surface  model,  three  struc¬ 
tures  are  modeled.  The  first  is  an  abstract  geometry  having  the  general 
shape  of  a  complex  jet  inlet  duct.  The  second  is  the  front  end  of  a  par¬ 
tial  jet  turbine.  And  the  final  is  a  model  of  a  missile.  Three  views  of  the 
first  model  are  given  in  figure  11.  The  second  model  is  given  in  figure  12. 
It  should  be  noted  that,  in  this  figure,  the  turbine  blades  are  not  simple 
plates  but,  rather,  possess  a  finite  thickness  away  from  each  edge.  A  full 
disk  of  36  blades  could  have  been  incorporated  but  would  have  required  a 
redimensioning  of  the  model-implementing  computer  code — more  on  this  in 
appendix  A.  The  final  model  is  given  in  figures  13  and  14.  To  demonstrate 
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the  role  of  the  topological  transformation  in  this  model,  the  first  figure  is  the 
missile  with  the  transformation  suppressed.  It  should  be  noted  that  all  tail 
surfaces  taper  from  an  edge  to  a  finite  central  thickness  and  that  the  wings 
have  a  complex  airfoil  shape  characteristic  of  subsonic  aircraft.  From  the 
first  two  models,  it  should  be  apparent  that  an  internal  inlet  duct,  turbine 
front  end,  exhaust  duct,  and  engine  outlet  could  have  been  incorporated  into 
this  model.  It  must  be  stressed  that  the  blades  of  figure  12  and  the  wing 
and  tail  surfaces  of  figures  13  and  14  are  not  modeled  separately  and  then 
attached  to  the  central  body  but,  rather,  are  designed  as  an  integral  part  of 
the  central  body.  This  is  the  only  reasonable  way  of  avoiding  the  inaccuracy 
inherent  in  the  attempted  mating  of  surfaces  with  complex  shapes. 

2.2.6  Concluding  Remarks 

The  equations  of  the  analytic  model  represent  the  best  perceived  compromise 
between  simplicity — with  all  its  benefits — and  the  need  to  model  complicated 
surfaces  in  the  most  direct  way  possible.  As  a  consequence,  the  most  produc¬ 
tive  upgrades  are  not  to  the  model  but  rather  to  the  software  implementation 
of  the  model.  For  instance,  one  useful  improvement  to  the  code  is  to  permit  a 
rotation  and  translation  in  three  dimensions  of  individual  model  independent 
structures.  Another  is  a  further  generalization  of  conditions  for  defining  the 
interval  over  which  the  topological  transformation  operates.  It  is  expected 
that  future  versions  of  the  code  will  contain  these  upgrades. 

The  key  to  the  proper  use  of  the  model  is  the  ability  to  define  an  appropriate 
set  of  cross  sections.  As  has  been  noted,  cross  sections  are  formed  by  joining 
limited  numbers  of  curve  segments.  Ideally,  one  would  hope  that  the  cross- 
section  function  could  define  the  sought-for  shapes  with  a  single  curve  and, 
in  lieu  of  this,  with  as  few  curves  as  possible.  To  construct  a  cross  section 
from  curve  segments,  it  is  very  helpful  to  define,  based  on  the  cross-section 
function,  a  set  of  curve  primitives.  One  obvious  choice  is  the  superellipse 
itself,  or  segments  thereof.  Figures  6  and  7  define  other  potential  candidates. 
Experience  indicates  that  a  potentially  fruitful  area  to  explore  for  defining 
curve  primitives  would  be  cross-section  functions  where  parameters  m  and  n 
are  not  set  equal  or  set  to  1.  Such  curves  do  not  have  the  simplicity  or  visual 
appeal  of  figure  7,  as  an  example,  but  do  appear  to  generate  potentially 
useful  curve  shapes. 
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Figure  1.  Superellipse. 
From  rectangle,  e  =  0.01, 
0.5,  1,  2,  4,  8. 


Figure  2.  Transition  func¬ 
tion  based  on  modified 
superellipse.  Curve  pa¬ 
rameters  are  listed  in 
appendix  B. 
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Figure  3.  Absolute  value 
of  first  derivative  of  tran¬ 
sition  function. 


Figure  4.  Absolute  value 
of  second  derivative  of 
transition  function. 
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Figure  5.  Curves  of  f(v)  =  given  in  legend. 
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m  =  n  =  1  uy  =  0  Uj  =  45 
€  =  0.01,0.5,1,2,4 


m=l  n  =  4  uy  =  u,  =  0 
£  =  0.01,0.5,1,2,4 


m=l  n  =  2  uv  =  u,  =  0 
€  =  0.01,0.5,1,2,4 


m  =  l  n  =  2  Uy  =  0  u,  =  45 
€  =  0.01,0.5,1,2,4 


Figure  6.  Examples  of  cross-section  function-step  1  equations. 
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m  =  n  =  1  uv-ux  =  u,  =  0  j'  =  2 
&i  =  l  6 j  =  3  63  =  4,8,16,32,64,128 


m  =  n=  l  tt„  =  u,  =  u7  =  0  j  =  2 
6j  =  l  62  =  7  63  =  4, 8, 16, 32, 64, 128 


m  =  n  =  1  tty  =  u,  =  0  u?  =  67.5  j  =  0.5 
61  =  1  62  =  7  63  =  200,400,800,1600,3200,6400 


m  =  n  =  1  uy  =  ut  =  u,  =  0  j  =  2 
61  =  1  62  =  -0.9  63  =  4, 8, 16, 32, 64, 128 


Figure  7a.  Examples  of  cross-section  function— step  2  equations. 
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m  =  n  =  1  uy  =  u,  =  uq  =  0  j  =  20 
&!  =  1  63  =  -0.95  6,  =  10 


m  =  n  =  1  uv  =  uM  =  u,  =  0  j  =  20 
61  =  2  62  =  —1  63  =  0.05 


m  =  n  =  1  uy  =  u,  =  0  u,  =  45  j  =  2  m  =  n  =  1  u„  =  u,  =  u,  =  0  j  =  20 

i>!  =  1  62  =  -0.9, 0.5, 1, 3, 7, 20  63  =  20  6X  =  10  62  =  -9  63  =  0.05 


Figure  7b.  Examples  of  cross-section  function-step  2  equations. 
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Figure  8.  Examples  of 
cross-section  function- 
step  4  equations.  Curve 
parameters  are  listed  in 
appendix  B. 


(d) 
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Figure  9.  Examples  of 
effect  of  size  function  equa¬ 
tions  on  right  circular 
cylinder  geometry.  Curve 
parameters  are  listed  in 
appendix  B. 


Figure  10.  Topologically 
transformed  rectangu¬ 
lar  cross-section  tube. 
Surface  parameters  are 
listed  in  appendix  B. 
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Figure  11.  Duct  model. 


Figure  12.  Turbine  front 
end  with  three  blades. 
Surface  parameters  are 
listed  in  appendix  B. 
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HU 


Figure  13.  Exploded  view 
of  missile  before  topo¬ 
logical  transformation 
(model  composed  of  five 
independent  structures). 


Figure  14.  Missile  af¬ 
ter  topological  transfor¬ 
mation.  Parameters  are 
listed  in  appendix  B 
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Appendix  A 

Source  Code  Listing 


Included  within  this  appendix  is  a  FORTRAN  source  code  listing  of  the 
analytic  surface  model  and  a  listing  of  a  code  for  implementing  only  the 
cross-section  function. 

If  the  user  wishes  to  increase  the  limits  on  the  various  parameters  of  the 
analytic  surface  model  code,  this  is  especially  easy  to  do.  The  code  limits 
on  parameters  such  as  number  of  curves  per  cross  section  or  number  of  cross 
sections  per  model  (as  defined  in  the  code  user’s  manual  section)  were  all 
selected  to  have  different  values.  This  allows  the  user  to  readily  determine 
the  dimension  statements  that  require  changes.  The  code  is  presently  di¬ 
mensioned  to  calculate  a  maximum  120  by  240  grid  of  surface  points  (120  in 
the  x  (longitudinal)  direction  and  240  transverse).  Again,  these  limits  are 
readily  changed  within  the  block  of  dimension  statements.  The  only  other 
potential  dimension  statement  upgrade  is  to  ngroup  (defined  in  the  user’s 
manual).  In  addition,  the  following  lines  of  the  surface  code  may  need  to  be 
changed  after  any  redimensioning: 

•  line  178  (do  212  jbj= 1,15) 

•  line  187  (do  202  jvv=id,2000,5) 

•  line  264  (do  110  ic=l,15) 

•  line  288  (do  112  ic=l,15) 

The  upper  limit  on  the  statements  of  items  1,  3,  and  4  reflects  the  maximum 
permissible  number  of  curve  segments  per  cross  section.  For  item  2,  the  limit 
of  2000  should  reflect  the  size  of  the  dimension  statement  for  array  ngroup. 
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c  ANALYTIC  SURFACE  MODEL  CODE 
c 

dimension  x(120) ,  xl(25) ,  yc(240,25),  zc(240,25) 
dimension  al(20),  a2(20),  a3(20) ,  a4(20),  a5(20) ,  nblock(lO) 
dimension  y(240,120).  z(240,120).  ns(10),  anx(10) 
dimension  cl(10),  c2(10),  c3(10) ,  c4(10),  c5(10) ,  xt (240 .120) 
dimension  an(15,25) ,am(15,25) ,alim(15,25) ,uy(15.25) ,uz(15,25) 
dimension  aj  (15,25)  ,bl(15,25)  ,b2(15.25)  ,b3{4r5,25)  ,a(15.25) 
dimension  b(15,25) ,psi(15,25) ,yt(15,25) ,zt(15.25) 
dimension  aa(15,25) ,bb(15.25) ,ang(15,25) .  npoint(lO) 
dimension  nseg(25) ,ytp(15,25) ,ztp(15,25) ,uq(15,25) 
dimension  czl(10),  cz2(10),  cz3(10) ,  cz4(10),  cz5(10) 
dimension  cyl(10),  cy2(10).  cy3(10),  cy4(10) ,  cy5(10) 
dimension  azl(20),  az2(20) ,  az3(20),  az4(20),  az5(20) 
dimension  ayl (20) ,  ay2(20),  ay3(20) ,  ay4(20),  ay5(20) 
dimension  ymin(15, 25) ,  ymax(15,25),  zmin(15,25),  zmax(15,25) 
dimension  npn( 15,25) ,  npnn(15,25),  yminn(15,2),  ymaxx(15,2) 
dimension  zminn( 15,2) ,  zmaxx(15,2),  ngroup(2000) 
dimension  h( 15,25) ,  hfl(15,25),  hf2(15.25),  hf3(15,25) 
dimension  hf4(15,25),  hf5(15,25),  hf6(15.25),  hf7(15,25) 
dimension  hf8(15,25),  hf9(15,25),  hbl(15,25),  hb2(15.25) 
dimension  hb3( 15,25) ,  hb4(15.25),  hb5(15,25),  hb6(15,25) 
dimension  hb7(15,25),  hb8(15,25),  hb9(15,25) 
dimension  hf 10(15,25) ,  hbl0(15,25),  ntop(10),  nwi(10) 
c 

namelist /lista/an.am.alim, uy ,uz ,aj ,bl ,b2 ,b3 ,a ,b,psi ,yt ,zt , 

&ang , nseg , npoint , aa , bb , y tp , ztp , uq . nblock , li st3 , 

&al , a2 ,a3 , a4 , a5 ,ns , anx.xl , c 1 , c2 ,c3 , c4 ,c5 .ncount ,mid , 

&czl , cz2 ,cz3 ,cz4 ,cz5 ,cyl ,cy2 ,cy3 ,cy4 ,cy5 , 

&azl ,az2,az3,az4,az5,ayl ,ay2,ay3,ay4,ay5, 

&h  ,hf  1  ,hf  2  ,hf  3  ,hf  4  ,hf  5  ,hf  6  ,hf  7  ,hf  8  ,hf  9  ,hf  10 , 

&hbl ,hb2 ,hb3 ,hb4 ,hb5 ,hb6 ,hb7 ,hb8 ,hb9 ,hblO .ngroup ,nwi ,ntop 
c 

c  treat  each  block  of  sections  defined  by  a  discontinuous  station  cut 
c  as  independent .. .ncount=  number  of  independent  blocks-  max. =5 
c  specify  all  station  cuts  for  each  independent  block 
c  ns=  number  of  model  3-D  segments  defined  by  station  cuts 
c  for  each  independent  block. . .max.  value=  20 
c  anx=  surface  grid  count  in  x  direction 
c 
c 
c 
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c 


read(3,lista) 
jv=l 
j  w=0 
kq=0 
kqq=0 
nsa=l 
88  i=0 

ia=l 

do  60  kz=l .ncount 
kzq=anx(kz) +1 . 1 
nsb=ns (kz) +nsa 
w=0. 

do  50  j=nsa,nsb-l 
w=xl ( j  + 1 ) -xl ( j ) + w 
50  continue 

w=w/ (anx(kz)-l . ) 
if (xl(nsa+l) .It .xl(nsa))  then 
signa=-l . 
else 

signa=+l . 
end  if 
jj=o 
jq=kq+l 

if(mid.ne.O)  then 
wa=xl(nsa)+w/2. 
go  to  2 
end  if 
wa=xl(nsa) 

2  do  62  j=l,kzq 

if (signa .gt .0. )  then 

if (wa.gt.xl(nsb))  go  to  63 

kq=l+kq 

x(kq)=wa 

wa=x(kq)+w*(l .-1 .e-05) 
go  to  62 
else 

if (wa.lt .xl(nsb))  go  to  63 

kq=l+kq 

x(kq)=wa 

wa=x(kq)+w*(l . -1 .e-05) 
end  if 
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62  continue 

63  continue 

c  cross-section  definitions: 

c  store  y  and  z  coordinates  of  surface  points  in  arrays  y(n,m),  z(n,m) 
c  where  m  defines  the  x  interval 
c  nseg=  number  of  curve  segments-  max. =15 
c  npoint=  number  of  plot  points 

c  ang=  angular  extent  of  each  segment  to  be  used  to  define  point 
c  distribution  between  segments  summation (ang) =360 
c  alim=  angular  stop  points  for  each  segment  for  the  sin  and  cos 
c  arguments  independent  of  parameters  am  and  an 
c  allow  for  axial  scaling  before  and  after  rotation 
c  each  segment  to  be  rotated  and  translated  independently 
do  1  j=nsa,nsb 
k=0 
npp=0 

do  3  n=l,nseg(j) 

np=npoint (kz) *ang (n , j ) /360 . +0 . 001 
npn(n, j)=np 

c  np  and  npn  are  number  of  plot  points  in  each  segment 
c  always  make  uy=  uz 

c  insure  that  summation  of  np*nseg=  constant  for  all  model  sections 
if(np.le.l)  go  to  5 
uu= (alim (n , j ) -uy (n , j ) ) / (np- 1 ) 
c  uu=  angular  increment  size 
u=0 . 

ymin(n, j)=l ,e06 
ymax(n, j)=-l .e06 
zmin(n, j)=l ,e06 
zmaxCn, j)=-l ,e06 
do  4  m=  1 ,  np 
k=l+k 

q=abs (sind(aj (n, j ) *u+uq(n , j ) ) ) 
eps=bl (n, j )  +  b2(n, j)*q**b3(n, j) 
atest=sind(am(n, j ) *u  +  uy(n,j)) 
btest=cosd(an(n, j)*u  +  uz(n,j)) 
aaa=sign(l . .atest) 
bbb=sign(l. .btest) 

yy=aaa*a (n , j ) *abs (sind(am(n , j ) *u+uy (n , j ) ) ) **eps+ytp (n . j ) 
zz=bbb*b (n . j ) *abs (cosd(an(n , j ) *u+uz (n , j ) ) ) **eps+ztp (n , j ) 
yc (k , j ) =aa(n , j ) * (yy*cosd(psi (n , j ) ) +zz*sind (psi (n . j ) ) ) +yt (n , j ) 
zc(k, j)=bb(n, j)*(yy*-sind(psi(n. j ) ) +zz*cosd(psi (n, j)))+zt(n, j) 
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if (yc(k, j) .It .ymin(n, j))  ymin(n, j )=yc (k, j ) 
if (yc(k, j) .gt.ymax(n, j))  ymax(n, j)=yc(k, j) 
if (zc(k, j) .It ,zmin(n, j))  zmin(n, j)=zc(k, j) 
if (zc(k, j) . gt.zmax(n, j))  zmax(n, j)=zc(k, j) 
c  y-zmin,  y-zmax  define  minimum  and  maximum  values  of  y  and  z 
c  for  each  segment  of  each  defined  cross  section .. .for  use 
c  with  topological  transform 
u=u+uu 

4  continue 

5  continue 
npp=np+npp 

3  continue 

write (* ,101)npp 

101  format (i9) 

1  continue 

do  64  j=nsa,nsb-l 
i=l  +  i 
c 

id=ia 

loop=0 

jb=j 

mb=l 

ib=0 

if (ngroup(ia) .eq.O)  go  to  100 

102  if (ngroup(ia) .ne . jb)  go  to  103 
ib=l+ib 

yminn(ib,mb)=l .e06 
ymaxx(ib,mb)=-l .e06 
zminn(ib,mb)=l .e06 
zmaxx(ib,mb)=-l .e06 
do  104  ma=ngroup(ia+l) ,ngroup(ia+3) 

if (ymin(ma, jb) .lt.yminn(ib.mb))  yminn(ib,mb)=ymin(ma, jb) 
if (ymax(ma , jb) . gt . ymaxx(ib.mb) )  ymaxx(ib,mb)=ymax(ma, jb) 
if (zmin(ma , jb) . It .zminn(ib ,mb) )  zminn(ib ,mb)=zmin(ma , jb) 
if (zmax(ma , jb) .gt .zmaxx(ib.rnb) )  zmaxx(ib ,mb)=zmax(ma , jb) 
npnn ( ib , j  b ) =npn (ma , j  b) +npnn ( ib , j  b) 
write ( * , 220) ib , j  b , npnn ( ib , j  b) , mb , 

&ia .ngroup(ia) .loop 
220  format (7i5) 

104  continue 

c  npnn  is  number  of  cross-section  points  for  each  cluster 
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ia=ia+5 

if (ngroup(ia) .eq. jb)  go  to  102 

if (nwi (kz) . ne .0 . and . j . ne . (nsb-1) . and . loop . eq. 1)  ia=id 
103  if (nsb-1 . eq. j )  then 

if (loop.eq. 1)  go  to  100 

201  ib=0 
mb =2 
jb=nsb 
loop=l 

do  212  j  b j  =  1,15 
npnn(jbj , jb)=0 
212  continue 

go  to  102 
end  if 

if (nwi(kz) .eq.O)  go  to  100 
if (loop.eq. 1)  go  to  100 
if(ib.eq.O)  go  to  100 
id=ia 

do  202  jvv=id,2000,5 
if (ngroup(jvv) .eq.nsb)  then 
ia=j vv 
go  to  201 
end  if 

202  continue 
write (* ,203) 

203  format (6h  error) 

100  continue 

c  ngroup  defines  clustering  of  cross-section  segments  for  use  by 
c  topological  transform  -  example  input:  1 , 1 , 1 ,3.0, 1 .4 , 1 ,8 ,0,2 . 1 ,2 . 8 , 
c  use  0  to  separate  clusters  . . .  cluster  pairs  define  array  elements 
c  for  cross-section  segments  -  must  be  no  gaps  in  any  cited 
c  cross  section  or  cited  independent  block 
c 

c  size  function  can  be  specified  for  y  and  z  axis  independently 
c  through  parameters  cyl ,czl ,cy2 ,cz2 , . . .  or  for  both  y  and  z 
c  axes  through  cl,c2,... 

c  similarly,  for  the  transition  function  as  ayl ,azl ,ay2 , az2 
c  or  al,a2, .  .  . 

do  65  j j=j j j ,kq 

xx=(x( j  j ) -xl(nsa) )/ (xl(nsb) -xl (nsa) ) 
xa= (x ( j  j ) -xl ( j ) ) / (xl ( j  + 1 ) -xl ( j ) ) 

jjj=jj 
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if(xa.gt.l.)  go  to  64 
if(xx.gt.l.)  go  to  64 

ayz=abs(azl(i))+abs(az2(i))+abs(az3(i))+abs(az4(i))+abs(az5(i)) 
if (ayz . It . 1 . e-06)  then 
If lag=0 
else 

If lag- 1 

go  to  12 
end  if 

eps=2 . / (al (i)+a2(i) *xa**a3(i) +a4 (i) *xa**a5(i) ) 
go  to  13 

12  epsz=2./ (azl(i)+az2(i)*xa**az3(i)+az4(i)*xa**az5(i)) 
epsy=2 . / (ayl (i) +ay2(i) *xa**ay3(i) +ay4(i) *xa**ay5(i) ) 

13  cyz=abs (czl (kz) ) +abs (cz2(kz) ) +abs (cz3(kz) ) +abs (cz4 (kz) ) + 
&abs(cz5(kz)) 

if (cyz . It . 1 . e-06)  then 
kf lag=0 
else 

kf lag=l 
go  to  10 
end  if 

s=cl (kz) +c2(kz) *xx**c3(kz)+c4(kz) *xx**c5(kz) 
go  to  11 

10  sy=cyl (kz) +cy2(kz) *xx**cy3 (kz) +cy4 (kz) *xx**cy5 (kz) 
sz=czl (kz) +cz2 (kz) *xx**cz3 (kz) +cz4 (kz) *xx**cz5 (kz) 

11  if (If lag. eq. 1) then 

epy=epsy 

epz=epsz 

else 

epy=eps 
epz=eps 
end  if 

ty=(l .0-xa**epy) ** (1 . /epy) 
tz=(l  .0-xa**epz)**(l  ./epz) 
kqq=l+kqq 
do  66  n=l ,k 
if (kflag.eq. 1)  then 
syy=sy 
szz=sz 
else 
syy=s 
szz=s 
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end  if 

y (n,kqq)=syy* ( (1 . -ty) *yc (n, j  +  1) +ty*yc (n, j ) ) 
z (n , kqq)=szz* ( ( 1 . -tz) *zc (n , j  +  1) +tz*zc (n , j ) ) 
c  ntop(n)=l  means  that  a  topological  transform  is  to  be 
c  performed  upon  independent  block  n 

c  to  include  topological  transform  of  lofted  model  simply  add 
c  h  parameters  to  namelist  input: 
if (ntop(kz) .eq.O)  go  to  80 
gf=0. 
gb=0. 
nx=0 

do  110  ic=l , 15 
nx=npnn(ic , j )  +nx 
if(nx.ge.n)  go  to  111 

110  continue 

111  icc=ic 

c  icc  determines  the  operational  cross-section  cluster 
ymm=ymaxx(icc , 1) -yminn(icc ,1) 
zmm=zmaxx(icc , 1) -zminn(icc , 1) 
if (ymm.lt . 1 .e-06)  then 
yf=.5 
go  to  115 
end  if 

yf=abs((y(n,kqq)-yminn(icc ,1))/ (ymaxx(icc , 1) -yminn(icc , 1) ) ) 

115  if (zmm.lt . 1 .e-06)  then 

zf=  .5 
go  to  116 
end  if 

zf=abs((z(n,kqq)-zminn(icc ,l))/(zmaxx(icc ,l)-zminn(icc ,1))) 

116  epf=hf5(icc . j)+hf6(icc , j ) *yf **hf 7(icc , j)+hf8(icc, j ) *zf **hf 9(icc , j) 
gf =hf 1 (icc , j ) *abs(l . -abs (hf 2(icc , j ) -yf ) ** (2 . /epf ) ) ** (0. 5*epf ) 

&  +hf 3(icc ,  j) *abs (1 . -abs (hf 4(icc , j ) -zf ) ** (2 . /epf ) ) ** (0.5*epf ) 

&  +hfl0(icc,j) 
if(mb.ne.2)  go  to  114 
nx=0 

do  112  ic=l ,15 
nx=npnn ( i c , nsb ) +nx 
if(nx.ge.n)  go  to  113 

112  continue 

113  icc=ic 

ymm=ymaxx(icc ,2)-yminn(icc ,2) 
zmm=zmaxx(icc ,2)-zminn(icc ,2) 
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if (ymm. It . 1 . e-06)  then 
yb= .  5 
go  to  117 
end  if 

yb=abs ( (y (n ,kqq) -yminn(icc ,2) ) / (ymaxx(icc ,2) -yminn(icc ,2) ) ) 

117  if (zmm. It . 1 . e-06)  then 

zb=  .5 
go  to  118 
end  if 

zb=abs ( (z (n ,kqq) -zminn(icc ,2) ) / (zmaxx(icc ,2) -zminn(icc ,2) ) ) 

118  epb=hb5(icc ,nsb) +hb6(icc ,nsb) *yb**hb7 (ice ,nsb)+ 

&hb8(icc ,nsb) *zb**hb9(icc ,nsb) 

gb=hbl (icc ,nsb) *abs (1 . -abs (hb2 (icc ,nsb) -yb) ** (2 . /epb) ) ** (0 ,5*epb) 
k  +hb3(icc ,nsb) *abs (1 . -abs (hb4 (icc ,nsb) -zb) ** (2 . /epb) ) ** (0 . 5*epb) 
k  +hbl0(icc ,nsb) 

114  continue 

c  limits  on  topological  transform  can  span  either  width  of  independent 
c  block  (specify  nwi=l)  or  width  of  segment  -  for  width  of  independent 
c  block:  repeat  hf  parameters  for  each  incorporated  cross-section, 
c  Need  include  hb  parameters  for  cross-section  nsb  only  for  nwi=l  or  0. 
if (nwi (kz) .eq. 1)  then 
xax=xx 
else 
xax=xa 
end  if 

xtt=gf *(1 . -xax**h(icc , j))+gb*xax**h(icc , j) 

gf=0. 

gb=0 . 

80  xt(n,kqq)=x(j j)+xtt 
xtt=0. 

c  nblock  allows  independent  blocks  to  be  stored  on  separate  groups  of 
c  fortran  logical  units-  nblock  defines  last  block  of  each  group 
if (kz . le .nblock( jv) )  go  to  70 
jv=l+jv 
j ww=j ww+3 
70  jx=60+jww 

if (list3 .ne . O.and .ngroup(l) . eq .0)  then 
write( jx, 192)x(kqq) ,y(n,kqq) ,z(n,kqq) 

192  format (3el5 .6) 
go  to  66 

else  if (list3 .ne .O.and.ngroup(l) .ne .0)  then 
write( jx, 192)xt (n.kqq) ,y(n,kqq) ,z(n,kqq) 
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go  to  66 
end  if 

if (ngroup(l) .eq.O.)  then 
wri te ( j  x , 190) x (kqq) 
write( jx+1 , 190)z(n,kqq) 
write (jx+2, 190) y(n, kqq) 
190  format (el6. 6) 
go  to  66 
else 

wr i t e ( j  x , 190) xt (n . kqq) 
write ( jx+1, 190) z(n, kqq) 
write ( jx+2, 190) y(n, kqq) 
end  if 

66  continue 
65  continue 
64  continue 
nsa=nsb+l 
60  continue 
stop 
end 
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c  CROSS-SECTION  CODE 
c 

dimension  an(90,5),  am(90,5),  alim(90,5),  uy(90,5),  uz(90,5) 
dimension  a j (90.5),  bl(90,5),  b2(90,5),  b3(90,5),  a(90,5) 
dimension  b(90,5),  psi(90,5),  yt(90,5),  zt(90,5),  y(720,5) 
dimension  z(720,5),  aa(90,5),  bb(90,5),  nbreak(90,5) ,  ang(90,5) 
dimension  nseg(5) ,  ytp(90,5),  ztp(90,5),  uq(90,5) 
name list /list a/an , am.alim.uy ,uz ,aj ,bl ,b2 ,b3 ,a ,b,psi ,yt ,zt , 

&ang ,nseg ,nplot .npoint ,aa ,bb .nbreak ,ytp ,ztp ,uq 
read(25 , lista) 

c  nseg=  number  of  curve  segments 
c  npoint=  number  of  plot  points 

c  ang=  angular  extent  of  each  segment  to  be  used  to  define  point 
c  distribution  between  segments  summation (ang) =360 
c  alim=  angular  stop  points  for  each  segment  for  the  sin  and  cos 
c  arguments  independent  of  parameters  am  and  an 
c  allow  for  axial  scaling  before  and  after  rotation 
c  each  segment  to  be  rotated  and  translated  independently 
do  1  j=l,nplot 
k=0 

j  3=1 

do  3  n=l,nseg(j) 
np=npo i nt  *  ang ( n , j ) / 360 . 

c  np  is  number  of  plot  points  in  each  segment 
c  always  make  uy=  uz 

uu= (alim (n , j ) -uy (n , j ) ) / (np- 1) 
c  uu=  angular  increment  size 
u=0. 

do  4  m= 1 , np 
k=l+k 

q=abs(sind(aj (n, j)*u+uq(n, j))) 
eps=bl(n,j)  +  b2(n, j)*q**b3(n, j) 
atest=sind(am(n, j)*u  +  uy(n,j)) 
btest=cosd(an(n, j)*u  +  uz(n,j)) 
aaa=sign(l . .atest) 
bbb=sign(l . .btest) 

yy=aaa*a (n , j ) *abs (sind (am (n . j ) *u+uy (n , j ) ) ) *  *eps+ytp (n , j ) 
zz=bbb*b (n , j ) *abs (cosd(an(n , j ) *u+uz (n , j ) ) ) **eps+ztp(n , j ) 
y(k, j)=aa(n, j)*(yy*cosd(psi(n, j))+zz*sind(psi(n, j)))+yt(n, j) 
z(k, j)=bb(n, j)*(yy*-sind(psi(n, j))+zz*cosd(psi(n, j)))+zt(n, j) 
c 

c  nbreak  used  to  set  number  of  segments  to  be  joined  to  form  individual 
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c  curves...  nbreak  defines  last  segment  in  each  curve 

c  j  term  allow  routine  to  be  used  for  multiple  cross  section  problem 

c 

if (n.le. nbreak (j j , j))  go  to  5 

5  kk=79+ j  j 

write (kk,6)z (k , j ) ,y(k, j) 

6  format(2el6.6) 
u=u+uu 

4  continue 
3  continue 
1  continue 
stop 
end 
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Appendix  B 

Source  Code  Input  Modules 


The  curve  parameters  for  figures  2  through  4  (as  defined  in  sect.  2.1.22 
in  main  body  of  text)  are  included  within  this  appendix.  In  addition,  the 
appendix  contains  a  listing  of  the  code  input  (namelist)  modules  for  a  num¬ 
ber  of  the  figures  of  this  report.  These  represent  sufficient  information  to 
demonstrate  the  performance  of  the  surface  code. 
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DATA  FOR  FIGURES  2  THROUGH  4 


CURVE 

al 

a2 

a3 

a4 

a5 

1 

0.01 

20.0 

0.5 

0.0 

— 

2 

0.01 

4.0 

0.3 

1.0 

4.0 

3 

0.01 

1.0 

4.0 

4.0 

1.0 

4 

0.01 

1.0 

10.0 

4.0 

4.0 

5 

0.01 

20.0 

50.0 

0.0 

— 

INPUT 

DATA 

SET  FIGURE 

00 

&li sta  npoint=360 , nseg=4 , ang=90 ,90,90,90, alim= 180 , 360 , 180 , 360 , 

uy=0, 180,0, 180, uz=0, 180,0, 180, uq=0, 0, 0, 0, bl=l, 1, 1, 1 ,b2=0, 0,0,0, 

b3=l ,1,1,1 ,am=3 ,1,3,1 ,an=l ,l,l,l,aj  =  l,l,l,l,a=.l,l,  .1,1, 

b= 1 , 1 . 1 , 1 , aa= 1 , 1 , 1 , 1 . bb= 1 , 1 , 1 , 1 , psi=90 . 90 . -90 , -90 , 

zt=l , 1 , -1 , -1 ,yt=0 ,0 ,0 ,0 ,ztp=0 ,0 ,0,0 ,ytp=0 ,0,0,0,nbreak=2 ,4 ,nplot=l , 

&end 


INPUT  DATA  SET  FIGURE  8B 

iklista  npoint=360 ,nseg=2 , ang=180 , 180 , alim=360 ,61 . 36 , uy=45 , 360 , 
uz=45,360,uq=0,0,0,0,bl=l , 1 ,b2=0,0,b3=l , 1 ,am=l , -1 ,an=l , -1 , 
aj=l ,1 ,a=l , .75 ,b=l , .75 ,aa=l , 1 ,bb=l , 1 ,psi=67 .5 ,59.32 ,zt=0,0, 
yt=0, .2789,ztp=0,0,ytp=0,0,nbreak=2,nplot=l , 

&end 


INPUT  DATA  SET  FIGURE  8C 

klista  npoint=240,nbreak=8 ,nplot=l ,nseg=8 , 
ang=30 ,60,30,60,30,60,30,60, 

alim=87 .45,180,177.45,270, 267 . 45 , 360 , 357 .45,90, 

uy=2 .55,-90,92.55,0, 182 . 55,90,272 .55,-180, 

uz=2. 55. -90, 92. 55, 0,182. 55, 90. 272. 55, -180, 

bl  =  l,2,l,2,l,2,l  , 2,b2=8*0,b3=8*l  ,8011=8*1  ,an=8*l  ,aj=8*l , 

bb=l,. 0106,1, .272,1,. 0106,1, ,272,aa=l, . 272 , 1 . .0106 , 1 , .272,1, .0106, 

a=. 169,1, .169,1, .169,1, .169,1, b=. 169,1, .169,1, .169,1, .169,1, 

psi=0,45 ,0 ,45 ,0 ,45 ,0,45 ,zt=0,0,0 , - .362,0,0,0, .362, 

yt=0, .362, 0,0,0, -.362, 0,0, 

icend 
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INPUT  DATA  SET  FIGURE  8D 


itlista  aj =1,1, 1,-1, 1,1,1, -1.1,1, 1 .-1 ,aa=90*l ,bb=90*l . 

am=l. -1.1,1 .1,-1 .1.1. 1,-1 tl.l,an=l,-l, 1,1, 1,-1, 1,1,1, -1,1,1, 

ang= 110,40,0,40,5,40,0,40,5,40,0,40, 

uy=27. 15, 90, 0.0, 7. 15, 90, 0,0, 17. 15. 90. 0,0, 

uz=27. 15, 90, 0,0, 7. 15, 90, 0,0, 17. 15. 90, 0,0, 

uq=0, 0,0, 90, 0,0, 0,90, 0,0, 0,90, 

alim=366 ,0,90,90,16,0,90,90,26,0,90,90, 

bl=l. 1.8, 2, 1.8, 1,1. 8. 2, 1.8, 1,1. 8, 2, 1.8, 

b2=0, .4,0, .4,0, .4.0, .4,0, .4,0, .4,b3=15*l, 

a=l ,2, .0141,2,1,2, .0141,2,1,2, .0141, 2, nseg=12 ,nplot=l . 

b= 1 , 2 , .0141,2,1,2, .0141,2,1,2, .0141 ,2 ,npoint=360, 

psi=0, 51, -39, 51, 0,61, -29, 61, 0,71, -19, 71 ,nbreak=12, 

yt=0, -1.154, .409, -1.134,0, -.694, 1.062, -.675,0, -.213, 1.683, -.195, 

zt=0 ,2.549,3.796,2.547,0,2.711,3.668,2.705.0.2.790,3.428.2.788. 

&end 


INPUT  DATA  SET 
FIGURE  9 A 

felista  ncount=l ,xl=0, 10,ns=l ,anx=60,npoint=30,nseg=l , 1 , 
ang=30*360,alim=30*360,bl=30*l  ,b3=30*l  ,am=30*l  ,an=30*l , 
a=30*l ,b=30*l ,aa=30*l ,bb=30*l ,nblock=l ,al=l ,a3=l ,a5=l , 
yt=30*1000,zt=30*0,czl=l ,cz3=l ,cz5=l , 
cyl=l ,cy2= .01 ,cy3= . 5 ,cy4=- .01 ,cy5=3 , 
itend 


INPUT  DATA  SET 
FIGURE  9B 

fclista  ncount=l ,xl=0, 10,ns=l ,anx=60,npoint=30,nseg=l , 1 , 
ang=30*360, 311111=30*360,61=30*1  ,b3=30*l  ,am=30*l  ,an=30*l , 
a=30*l ,b=30*l ,aa=30*l ,bb=30*l ,nblock=l ,al=l ,a3=l ,a5=l , 
yt=30*1000,zt=30*0 ,czl=l ,cz2=2 ,cz3=2 ,cz4=0,cz5=l , 
cyl=l ,cy2=.01,cy3=.5,cy4=-.01 ,cy5=3, 

&end 
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INPUT  DATA  SET 
FIGURE  9C 

tlista  ncount=l ,xl=0, 10 ,ns=l ,anx=60,npoint=30,nseg=l , 1 , 
ang=30*360 , alim=30*360 ,bl=30*l ,b3=30*l ,am=30*l ,an=30*l , 
a=30*l ,b=30*l ,aa=30*l ,bb=30*l ,nblock=l , al=l , a3=l ,a5=l , 
y t =30*0 ,zt =30*0 , czl=4 , cz2=-3 , cz3= . 5 , cz4=0 , cz5=l , 
cyl=l ,cy2=3,cy3=2,cy4=0,cy5=l , 
tend 


INPUT  DATA  SET  FIGURE  10 

tlista  ncount=2,nblock=l ,2 ,xl=2 ,4 ,4 ,6 ,ns=l , 1 ,anx=50,50,npoint=120, 120, 

ang=4*90, 11*0, 90, 45. 45, 90, 45, 45, 9*0, 90. 45. 45, 90. 45. 45. 9*0, 4*90, 

alim=90 ,180.270,360,11*0,90,135.180,270.315,360,9*0. 

90,135.180,270,315,360,9*0,90,180,270,360.11*0, 

uy=0, 90, 180, 270, 11*0, 0,90, 135, 180, 270, 315, 9*0, 

0,90.135,180,270,315,9*0.90,180,270.11*0. 

uz=0, 90. 180, 270, 11*0, 0,90, 135, 180, 270, 315, 9*0, 

0.90,135,180,270,315,9*0,90,180,270,11*0. 

nseg=4 ,6 , 6 ,4 ,bl=60*2 ,b3=60*l ,am=60*l ,an=60*l ,aj=60*l , 

a=60* 1 ,b=60* 1 , aa=60*l ,bb=60* 1 . al=l , 1 , a3=l , 1 , a5=l , 1 , 

cl=l,l,c3=l,l,c5=l,l ,list3=0 ,psi=60*45 ,h=60*l , 

hfl=30*-l,0,l,1.0,l,l,hf 2=30*0, 0,0, 1.0. 1,0, 

hf 5=30*1 . 1 , 1 , 1 , 1 . 1 , 1 ,hf6=30*2 ,hf 7=30*1 , 1 , 1 , 1 . 1 . 1 , 1 . 

hf 9=30*1, l.l,l,l,l,l,hbl=0,- 1,-1, 0,-1, -1,9*0, 15*0, 15*0, 

hb2=0.0,l,0.1.0.hb5=l. 1,1, 1,1, 1,9*0, 15*0, 15*2, 

hb7=l , 1 , 1 , 1 . 1 , 1 ,9*0, 15*0, 15*1 ,hb9=l , 1 . 1 , 1 , 1 , 1 .9*0 , 15*0 . 15*1 . 

ngroup= 1,1, 1,4, 0,2, 1,2, 1,0, 2, 2, 2, 2, 0,2, 3, 2, 3, 0,2, 4, 2, 4,0, 

2, 5, 2, 5, 0.2, 6, 2, 6, 0.3, 1.3, 1,0, 3. 2. 3, 2. 0.3, 3, 3. 3, 0.3. 4, 3, 4,0, 

3, 5, 3, 5, 0,3, 6, 3, 6, 0,4,1 ,4,4,0,nwi=l ,ntop=l ,1 , 
tend 

INPUT  DATA  SET  FIGURE  12 

tlista 
ncount=2 , 

xl=0 .1.732,1.732,1.882.2.032, 

ns=l ,2 , 

anx=20,20. 
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npoint=240 , 240 , 
nseg=l ,1,9,12,9, 

ang=360, 14*0, 360, 14*0, 78. ,45,45,6,45,45,6,45,45,6*0, 

78. . 45.0.45.6.45.0.45.6.45.0.45.3*0, 

78. . 45.45.6.45.45.6.45.45.6*0. 
uy=15*0, 15*0, 20, 90, 0,0, 90, 0,10, 90, 0,6*0, 
27.15,90,0.0,7.15,90,0,0,17.15,90,0,0,3*0, 
32,90,0,12,90,0,22,90,0,6*0, 

uz=15*0, 15*0, 20, 90, 0,0, 90, 0,10, 90, 0,6*0, 

27.15,90,0,0,7.15,90,0,0,17.15,90,0.0.3*0, 

32.90,0,12,90,0,22,90,0,6*0, 

alim=360 , 14*0 , 360 , 14*0 , 360 ,0,90,10,0,90,20,0.90 ,6*0 , 

366,0,90,90,16,0,90,90,26,0,90,90,3*0, 

372.0,90,22,0,90,32,0,90,6*0, 

bl=l  ,14*0,1, 14*0 ,1,1. 8, 1.8, 1,1. 8, 1.8, 1,1. 8.1. 8,6*0 , 

1,1. 8, 2, 1.8, 1,1. 8, 2, 1.8, 1,1. 8, 2, 1.8, 3*0, 

1,1. 8, 1.8, 1,1. 8, 1.8, 1,1. 8, 1.8, 6*0, 
b2=15*0, 15*0,0, .4, .4,0, .4, .4,0, .4, .4,6*0, 

0,. 4,0,. 4,0.. 4,0,. 4,0,. 4,0, .4,3*0, 

0,. 4, .4,0, .4, .4,0, .4, .4,6*0, 
b3=75*l , 

am=30* 1,1, -1,1,1, -1.1. 1,-1. 1.6*1, 

1,-1, 1.1. 1,-1, 1,1, 1,-1, 1,1, 3*1, 

1,-1, 1,1, -1,1, 1,-1, 1,6*1, 
an=30*l, 1,-1, 1,1, -1,1, 1,-1. 1.6*1, 

1,-1, 1.1, 1,-1, 1.1. 1,-1, 1.1, 3*1, 

1,-1, 1,1, -1,1, 1,-1, 1.6*1, 
a j =30* 1,1, 1,-1, 1.1, -1,1. 1,-1. 6*0 , 

1.1. 1, -1, 1,1, 1,-1, 1,1,1. -1,3*0, 

1.1,  -1,1, 1,-1. 1,1. -1.6*0, 

a=l, 14*0, 1,14*0, 1,2, 2, 1,2, 2, 1,2, 2. 6*0, 

1.2. . 0141.2.1.2. .0141. 2. 1.2.. 0141. 2. 3*0, 

1 ,2,2,1 ,2,2,1 ,2, 2, 6*0, 

b=l, 14*0, 1,14*0, 1.2. 2, 1.2, 2. 1,2, 2, 6*0, 

1.2. . 0141.2.1.2. .0141.2.1.2. .0141.2.3*0. 

1,2, 2, 1,2, 2, 1,2, 2, 6*0, 

aa=75*l , 
bb=75*l . 

psi=15*0, 15*0, 0,45, 45, 0,55, 55, 0,65, 65, 6*0, 
0,51,-39,51.0.61,-29.61,0.71.-19,71.3*0, 

0,57,57,0,67,67,0,77,77, 

yt=15*0, 15*0, 0,-1. 414, -1.414,0, -.9734, -.9734.0, -.5032. -.5032. 6*0, 
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0,-1.154, .409, -1.134,0, -.694 ,1.062, -.675,0 ,-.213, 1.683, -.195, 3*0, 
0,-. 881, -.881,0, -.406, -.406,0, .080, .080, 

zt=15*0, 15*0,0,2.414,2.414,0,2.623,2.623,0,2. 752,2.752,6*0 , 

0,2. 549 ,3.796,2. 547 ,0,2.711,3. 668 ,2.705,0,2.790,3. 428 ,2.788,3*0, 
0.2.655,2.655,0.2.768,2.768,0,2.797,2.797,0, 
uq=15*0, 15*0, 0,0, 90, 0,0, 90, 0.0, 90, 6*0, 
0,0,0,90,0,0,0,90,0,0,0,90,3*0, 

0,0,90,0,0,90,0,0,90,6*0, 

nblock=l ,2 , 

list3=0, 

al=2,2,2, 

a2=0,0,0, 

a3=l ,1,1, 

a4=0,0,0, 

a5=l  ,1,1, 

cl-0.1, 

c2=l ,0, 

cS-l.l. 

c4=0,0, 

c5=l , 1 , 

&end 


INPUT  DATA  SET  FIGURE  14 
(fclista 

ncount=5 ,nblock=l ,2 ,3 ,4 ,5 ,ns=2 ,2 , 1 ,3,2, 

xl=l. 5, 2, 3. 9, 3. 9, 4. 07. 4. 36, 4. 36, 5. 46, 5. 46. 5. 96, 6. 33, 6. 86, 6. 86, 7. 1,7. 27. 

anx=42 ,8 ,19,24,7 ,npoint=60, 180,60,120,180, 

nseg=l, 1,1, 4, 6, 4, 1,1, 8, 8, 8. 8, 12, 8. 12, 

ang=360, 14*0,360, 14*0,360, 14*0, 

60,120,60,120,11*0.6*60.9*0,60,120,60,120,11*0, 

360,14*0,360,14*0, 

240,6,18,18.36,18,18,6,7*0,240,6,18,18,36.18,18,6,7*0, 
240,6,18,18.36,18,18,6.7*0.240,6,18,18,36.18,18,6,7*0. 
12*30,3*0,30,60,30,60,30,60.30,60,7*0,12*30,3*0, 
alim=360 , 14*0 , 360 . 14*0 , 360 , 14*0 , 

179,270,361.90.11*0,173.5,90.270.360.360,90,9*0,180,270,360,90,11*0. 

360,14*0,360,14*0, 

270,180,180,180,304.5,360,90,360,7*0,235.53,180,180,180,304.5.360,90. 

360.7*0,235.53.235.6,235.7,235.8,304.1,304.2,304.3.304.47.7*0, 

235 .53,235.6,235.7,235.8.304.1,304.2,304.3, 304 . 47 , 7*0 . 


53 


90,360,270 ,180,90,0,270, 180,90,360,270, 180,3*0,87.45 ,180, 177.45,270, 
267 .45,360,357.45,90,7*0,90,360,270, 180,90,0,270, 180,90,360,270, 
180,3*0, 

uy=15*0, 15*0 ,15*0, 

1,0,179,-180,11*0,6.5,0,90,180,180.0,9*0.0.0,180,-180.11*0. 

0,14*0,0,14*0, 

-90,90,90,90,235.5,270,0,270,7*0,-55.53,90.90,90,235.5.270.0,270,7*0. 

-55.53,235.53,235.6,235.7,235.8,304.1,304.2,304.3,7*0, 

-55.53,235.53,235.6,235.7,235.8,304.1,304.2,304.3,7*0. 

0,270,360,90,0,90,180,90,180,270,180,270,3*0,2.55,-90,92.55,0, 

182.55.90.272.55, -180,7*0,0,270,360,90,0.90.180,90.180,270,180, 
270,3*0, 

uz=15*0, 15*0, 15*0, 

1,0,179,-180,11*0,6.5,0,90,180,180,0,9*0,0,0,180,-180,11*0. 

0,14*0,0,14*0, 

-90,90.90,90,235.5.270.0,270,7*0,-55.53,90.90.90.235.5,270,0.270.7*0. 

-55.53,235.53,235.6,235.7,235.8.304.1,304.2,304.3,7*0, 

-55.53,235.53,235.6,235.7,235.8,304.1.304.2.304.3,7*0, 

0,270,360,90,0,90,180,90,180,270,180,270,3*0,2.55,-90,92.55.0,182.55, 

90.272.55, -180,7*0,0,270,360,90,0,90.180,90,180,270,180.270,3*0, 
bl=l ,14*0,1 ,14*0,1 ,14*0, 

1,2,1,2,11*0,1,2,2,1,2,2, 9*0 ,1,2,1,2,11*0, 


1,14*0,1,14*0, 

1,2, 1.2, 1,2, 1,2, 7*0, 1,2, 1,2, 1,2, 1,2, 7*0. 8*1, 7*0, 8*1, 7*0, 

1,2, 2, 1,2, 2, 1,2, 2, 1,2, 2. 3*0, 1,2, 1,2, 1,2, 1.2, 7*0, 1,2, 2, 1,2, 2, 1,2, 2, 
1.2, 2, 3*0, 


b2=15*0, 15*0, 15*0, 

15*0, 15*0, 15*0, 

15*0,15*0, 

15*0,15*0 , 15*0,15*0, 

15*0 , 15*0 , 15*0, 

b3=l, 14*0, 1,14*0, 1,14*0, 

1,1,1.1,11*0.6*1,9*0,1,1,1,1,11*0, 

1,14*0,1,14*0, 

8*1 ,7*0, 8*1 ,7*0, 8*1 ,7*0, 8*1 .7*0, 
12*1,3*0,8*1,7*0,12*1.3*0, 


am=l , 14*0 , 1 , 14*0 , 1 , 14*0 , 
1,1,1,1,11*0,6*1,9*0,1,1,1,1,11*0, 

1,14*0,1,14*0, 

1,1. 3. 1.1. 1,3, 1.7*0. 1.1. 3, 1.1, 1.3. 1.7*0, 8*1, 7*0, 8*1, 7*0, 

12*1,3*0,8*1,7*0,12*1,3*0, 

an=l ,14*0,1, 14*0, 1,14*0, 
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1,1,1,1,11*0,6*1 ,9*0, 1,1,1,1,11*0, 

1.14*0,1,14*0, 

8*1 ,7*0, 8*1 ,7*0, 8*1 ,7*0, 8*1 ,7*0, 

12*1,3*0,8*1,7*0,12*1,3*0, 
aj=l , 14*0,1 , 14*0,1 , 14*0, 

1,1.1,1,11*0.6*1.9*0.1.1,1,1,11*0, 

1,14*0,1,14*0, 

8*1 ,7*0, 8*1 ,7*0, 8*1 ,7*0, 8*1 ,7*0, 

12*1,3*0,8*1,7*0,12*1,3*0, 
a=l, 14*0, 1,14*0, 1,14*0, 

.265,1, .265,1,11*0.1. .735,1,1,1, .735,9*0,1,1,1,1.11*0, 

.265,14*0, .265,14*0, 

.265, .0354, -.015, .0707, .265, .0707, .015, .0354,7*0, .265, .0001 ,- .0001 , 
.0707, .265, .0707, .0001 , .0001 , 7*0 , 8* . 265 , 7*0 ,8* . 217 , 7*0, 

.217, .238, .238, .217, .238, .238, .217, .238, .238, .217, .238. .238,3*0, 
.169,1, .169,1, .169,1,. 169. 1,7*0, .12, .307, .307, .12, .307, .307, 

.12. .307, .307, .12, . 307 , . 307 , 3*0 , 
b=l, 14*0, 1,14*0, 1,14*0, 

.265,1, .265,1,11*0,1, .735,1,1,1, .735,9*0,1,1,1,1,11*0, 

.265,14*0, .265,14*0, 

.265, .0354, .15. .0707, .265, .0707, . 15, .0354,7*0, .265, .0001 , .0001 , .0707, 
.265, .0707, .0001, .0001 ,7*0,8* . 265 ,7*0,8* . 217 ,7*0, 

.217, .238, .238, .217, . 238 , . 238 , . 217 , . 238 . . 238 . . 217 , . 238 . . 238 . 3*0 . 

.  169 , 1 , . 169 . 1 , . 169 , 1 , . 169 , 1 , 7*0 , . 12 , . 307 , . 307 , . 12 , . 307 , . 307 , . 12 , 

.307, .307, .12, .307, .307,3*0, 
aa=. 0001, 14*0, .265,14*0, .265,14*0, 

1. . 0001.1. .0001.11*0, .265.1. .0071, .265, .0071.1,9*0. . 265 . .0001 . . 265 . 
.0001,11*0, 

1,14*0,1,14*0, 

8*1, 7*0, 8*1 ,7*0, 8* 1,7*0, 8*1 ,7*0, 

12*1,3*0,1, .272,1. .0106,1, .272,1, .0106.7*0.12*1.3*0, 
bb= . 0001 , 14*0 , . 265 , 14*0 , . 265 , 14*0 , 

1. . 636.1. .636.11*0, .265,1, .735, .265, .735.1.9*0, .265, .735, .265, .735, 
11*0, 

1,14*0,1,14*0, 

8*1 ,7*0, 8*1 ,7*0, 8*1 ,7*0, 8*1 ,7*0, 

12*1.3*0,1, .0106,1. .272.1, .0106,1, .272,7*0,12*1,3*0. 
psi=15*0, 15*0, 15*0, 

0,45.0,45,11*0,0,46.1 ,45,0,45,43.9,9*0,0.45,0,45,11*0, 

0.14*0,0,14*0. 

0,45,0.45,0,45,0,45,7*0,0,45,0,45,0,45,0,45,7*0,8*0,7*0,8*0,7*0, 

0,45,45,0,45,45.0.45,45.0,45,45,3*0,0,45,0,45,0,45,0,45.7*0, 
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0,45,45,0,45,45,0,45,45,0,45,45,3*0, 
zt=15*0, 15*0, 15*0, 

0.-.715.0, .715, 11*0,0, -.775, -.785,0, .785, .775, 9*0,0, -.785,0, .785,11*0, 
0,14*0,0,14*0, 

0, .025,0, -.1,0, .1,0, -.025, 7*0,0, -.15, -.15, -.1,0, .1, .15, .15,7*0, 

8*0, 7*0, 8*0, 7*0, 

0,-. 1685, -.1685,0, -.3855, -.3855,0, .1685, .1685,0, .3855, .3855,3*0, 

0,0,0, -.362, 0,0,0, .362, 7*0,0, -.217,- .217,0,- .337, -.337,0, 

.217, .217,0, .337, .337,3*0, 
yt=15*0, 15*0, 15*0, 

0, .005,0, .005, 11*0,0, -.5. .005,0, .005, -.5, 9*0, 0,0, 0,0, 11*0, 
0,14*0,0,14*0, 

0 , - . 29 , - . 315 . - . 365 , - . 1965 , - . 365 , - . 315 . - . 29 , 7*0 .0 , - . 218 , - . 218 , - . 268 . 
-.1,-. 268, -.218, -.218, 7*0, 8*0, 7*0, 8*0, 7*0, 

0, .3885, .3885,0, -.1685, -.1685,0, -.3855, -.3855,0, .1685, .1686,3*0, 

0, .362, 0,0,0, -.362, 0,0, 7*0,0, .337, .337,0, -.217, -.217,0, 

-.337, -.337,0, .217, .217,3*0, 
al=3,2,4, .5.2. .0001 ,2, .001 ,2,2, 
a2=0 ,0,1, 2. 5, 0,2, 0,2, 0,0, 
a3=l, 1,1, 1,1, .2,1, .3,1,1, 

34=0,0,0,0,0,0,0,0.0,0, 

35=1,1,1,1,1,1,1,1,1,1, 
c 1=5*1 ,c3=5*l ,c5=5*l , 

hf 1=180*0,0, .18, .18,0,0,0,0, . 18, .18,0,0,0,3*0,0, .18,0,0,0, . 180,0, 

0, 7*0,0, .18, .18,0,0,0,0, .18, .18,0,0,0,3*0, 

hf 2=180*0, 0, 1 ,1,0, 0,0, 0,0, 0,0, 0,0, 3*0, 0.1, 0,0, 0,0, 0,0, 7*0, 0,1 , 

1,0, 0,0, 0,0, 0,0, 0,0, 3*0, 

hf 3=45*0,0, .12,0, .12,11*0,0, .12, .12,0, .12, .12,9*0, 

0, .12,0, .12,11*0,90*0,0,0,0.0. . 18 , . 18 ,0 ,0 ,0 ,0 , . 18 , . 18 , 3*0 . 

0,0,0,. 18, 0,0,0, .18,7*0,0,0,0,0, .18, . 18 ,0 .0 ,0 ,0 . . 18 , . 18 . 3*0 , 
hf 4=45*0, 0,0, 0,1 ,11*0,0,0,0,0,1 ,1 ,9*0, 0,0, 0,1 ,11*0,90*0, 

0,0, 0,0, 0,0, 0,0, 0,0, 1,1, 3*0, 0,0, 0,0, 0,0, 0,1, 7*0. 

0,0, 0,0, 0,0, 0,0, 0,0,1 ,1 ,3*0, 

hf 5=45*2, 2, 2, 2, 2, 11*2, 2, 2, 2, 2, 2, 2, 9*2, 2. 2, 2, 2, 11*2, 90*0, 
12*2,3*0,8*2,7*0,12*2,3*0, 
hf 7=225*1 , 
hf 9=225*1 , 

hb3=45*0 , 0 , - . 08 , 0 , - . 08 , 1 1 *0 , 0 , - . 08 , - . 08 , 0 , - . 08 , - . 08 , 9*0 , 

0,-. 08,0, -.08. 11*0, 

hb4=45*0,0,0,0,l ,11*0,0,0,0,0.1,1,9*0,0,0,0,1 ,11*0, 
hb5=45*2, 2, 2, 2, 2, 11*2, 2, 2. 2, 2. 2, 2, 9*2, 2, 2, 2, 2, 11*2, 90*0, 45*2, 
hb7=225*l , 
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hb9=225*l , 
h=225*l , 


ngroup=4 ,1,4, 1,0, 4, 2, 4, 2, 0.4, 3, 4, 3, 0,4, 4, 4, 4, 0,5, 1,5, 1,0,5, 2, 5, 2.0, 

5, 3, 5, 3, 0,5, 4, 5, 4, 0,5, 5, 5, 5, 0,5, 6, 5, 6, 0,6, 1,6, 1,0, 6, 2, 6, 2, 0,6, 3, 6, 3,0, 

6,4,6,4.0,13,1,13,1,0,13,2.13,2,0,13,3,13.3,0,13,4,13,4.0, 

13.5,13,5,0,13,6,13,6,0,13,7,13,7,0.13,8.13,8,0,13,9.13,9,0. 

13,10,13,10,0,13,11,13,11,0,13,12,13.12,0,14.1,14,1,0,14.2,14.2.0, 

14,3,14,3,0,14,4.14.4,0,14,5,14,5,0,14.6,14,6,0,14,7.14.7.0. 

14,8,14,8,0,15,1,15,1,0,15,2,15,2,0,15,3,15,3,0,15,4,15,4,0, 

15,5,15,5,0,15,6,15,6,0,15,7,15,7,0,15,8,15,8,0,15,9,15,9,0, 

15,10,15,10,0,15,11,15,11,0,15,12,15,12,0, 

nwi=0, 1 ,0,0,1 ,ntop=0, 1 ,0,0,1, 


&end 
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